Interacting Growth Walk - a model for hyperquenched homopolymer glass? 
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We show that the compact self avoiding walk configurations, kinetically generated by the recently 
introduced Interacting Growth Walk (IGW) model, can be considered as members of a canonical 
ensemble if they are assigned random values of energy. Such a mapping is necessary for studying the 
thermodynamic behaviour of this system. We have presented the specific heat data for the IGW, 
obtained from extensive simulations on a square lattice; we observe a broad hump in the specific 
heat above the ^-point, contrary to expectation. 
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Linear polymers in a poor solvent are known [1] to as- 
sume globular configurations below a tricritical tempera- 
ture Tg, called the 0-point. These globules acquire denser 
minimum energy configurations at lower temperatures. 
In the case of random heteropolymers, the 'quenched' 
random interactions between the constituent monomers 
frustrate the evolution of the globules towards their min- 
imum energy configurations. They arc thus forced to 
freeze into higher energy configurations (local minima). 
In fact, the heteropolymer globules serve as 'toy models' 
for protein folding phenomenon [2] . It has been shown re- 
cently [3] that even homopolymer globules can freeze into 
glassy states, due to a self-generated disorder brought 
about by the competing interactions and chain connec- 
tivity during the cooling process. In this sense, the freez- 
ing of a homopolymer globule is said to be analogous to 
that of a structural glass. 

In a Monte Carlo study of this freezing process, we 
may choose a configuration from a canonical ensemble of 
Interacting Self Avoiding Walks (ISAW) [4] which repre- 
sents a linear polymer in equilibrium with a thermal bath 
at a temperature T (say, > Tg). Then, using a standard 
dynamical algorithm [5] , we may relax the chosen config- 
uration at a temperature preset {i.e., quenched) to a de- 
sired value less than Tg; deeper the quench, more difficult 
and time consuming it would be to realize a globular con- 
figuration. On the other hand, the Interacting Growth 
Walk (IGW) [6] is a simpler but more efficient algorithm 
for generating compact or globular Self Avoiding Walks 
(SAW); they are generated, step by step, by sampHng 
the locally available sites with appropriate Boltzmann 
factors, ea;p(/3(3n^^eo), where (3^^ is the 'growth' tem- 
perature, n'^pj{l < m < z — 1) is the number of non- 
bonded nearest neighbour (nbNN) contacts the site m 
will make, if chosen, on a lattice of coordination number 
z and — eo is the attractive energy associated with any 
nbNN contact. 

In this paper, we show that these kinetically generated 
IGWs represent the frozen configurations of a homopoly- 



mer globule with a self-generated disorder. Contrary to 
expectation, our simulations on a square lattice indicate 
an excess specific heat, characterizing these frozen states, 
above the 0-point. In fact, this simple model demon- 
strates that a meaningful statistical mechanical descrip- 
tion of an irreversible growth process involves an element 
of self-generated disorder brought about by ergodicity- 
breaking of the system. 

The growth of an IGW starts by first " occupying" an 
arbitrarily chosen site Tq of a regular c?-dimensional lat- 
tice of coordination number z whose sites are initially 
"unoccupied" (by monomers). The first step of the walk 
is taken in one of the z available directions by choosing an 
"unoccupied" nearest neighbours (NN) of rg, say ri, at 
random and with equal probability. Let the walk be non- 
reversing so that it has a maximum of z — 1 directions to 
choose from for the next step. Let {r™ \ m — 1, 2, Zj} 
be the "unoccupied" NN's available for the jth step of 
the walk. If Zj =0, the walk cannot grow further because 
it is geometrically "trapped". It is, therefore, discarded 
and a fresh walk is started from Tq. If Zj ^ 0, the walk 
proceeds as follows: 

Let n^NNij) be the number of nbNN sites of r™. Then, 
the probability that this site is chosen for the jth step is 
given by. 
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where the summation is over all the zj available sites. 
At "infinite" temperature {Pg = 0), the local growth 
probability Pm(rj) is equal to 1/zj and thus, the walk 
generated will be the same as the Kinetic Growth Walk 
(KGW)[7]. However, at finite temperatures, the walk 
will prefer to step into a site with more nbNN contacts. 
We have illustrated this local growth rule in Fig. 1(a) for 
IGW on a square lattice. Lower the growth temperature, 
less is the attrition (see the inset of Fig. 2) that the walk 
suffers while also being able to grow into more compact 
configurations. Moreover, it has been shown [6] that a 
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6'-point for this walk exists, and that the walk belongs to 
the same universality class {i.e., has the same values of 
the universal exponents, v and 7) as the SAW above, at 
and below the ^-point. 
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FIG. 1. A simple illustration of the IGW algorithm for 
generating walks from the origin, denoted by the open cir- 
cle, at a given growth temperature, Pq^ ■ (a) The sites A, 
B and C are available for making the fifth step. Choos- 
ing the site A will lead to one nbNN contact, whereas 
choosing the sites B or C will lead to none. Hence, 
the sites A, B and C will be chosen with probabilities 
e''<3/(2 + e^G^), 1/(2 + e^^) and 1/(2 + e^^) respectively, 
(b) The probability of growing this configuration is given by 
Pb = (l/4)(l/3)^(l/2)2(e^°/[2 + e'5G])2(e'5«/[l + e^^]). (c) 
The probability of growing this configuration, which is iden- 
tical to (b), is given by pc = {l/4){l/3)^ {e^^'=' /[2 + e^'^^]). 

Wo have repeated the IGW simulations on a square 
lattice for walks upto N = 8000, much longer than re- 
ported in ref.[6] and with better statistics. In Fig. 2, we 
have shown the TV-dependence of the exponent, v{N), 
obtained from the mean squared radius of gyration data, 
for various values of /3g in the range 3 to 10. We have 
estimated the asymptotic values of this exponent as sim- 
ple polynomial extrapolations of these v{N) values, and 
presented them in Fig. 3, along with also those obtained 
for /3 = 0, 1, 1.5 and 2 from the earlier data reported in 
ref.[6]. 

The transition from the SAW phase {v = 3/4) to the 
collapsed walk phase {v = 1/2) seems to be taking place 
over a narrow range of (ia values (~ 3.5 < (3g <~ 5.0), 
but this could still be due to limitations of our numerical 
work. The asymptotic estimates of v could improve not 
only with longer walks but also with larger number of 
successful walks, and this could rcsidt in narrower tran- 
sition regime. The ^-point for the IGW corresponds to 
a growth temperature given by ~ 4.5, which is close 
to our earlier value (^ 4) [6]. Thus, we see that IGW 
has all the three distinct phases (extended, ^^-point and 
collapsed) of SAW, realizable by tuning the growth tem- 
perature I3q^. 

However, the IGW does not represent a homopolymer 
in equilibrium with its environment at some bath temper- 
ature. Because, the set of all A^-stcp IGWs generated at 
a given growth temperature, Zigw{N; Pg), is not equiv- 
alent to the canonical ensemble of ISAWs, Zisaw{N; /3), 



for some bath temperature l3~'^. For example, in Fig. 1(b) 
and 1(c), we have shown two identical configurations 
which are expected to occur with the same probability 
in a canonical ensemble, but arc in fact grown with dif- 
ferent probabilities. This is a consequence of the fact 
that the local growth probability, pj{vj), of making the 
jth step to a site Yj depends on all the previous sites 
visited. Hence, the probability of generating an IGW 
configuration, C = {ro, ri, r^, ...}, has to be written as 
r^; i"o, Tj-i). Nonetheless, 

there must be a correspondence between the kinetically 
generated IGW and the canonical ISAW, especially be- 
cause the former can be tuned to belong to the same 
universality classes as the latter. 
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FIG. 2. The trend towards the asymptotic values of the 
exponent, v, for various values of /3gs (=3.0, 3.5, 3.7, 3.8, 
3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 5.0 and 10.0, from top to 
bottom. Inset: Semi-logarithmic plot of the attrition con- 
stant as a function of /3g- The data seem to suggest a form, 
\iGw oc exp{—al3G), where a is a constant. 
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FIG. 3. The collapse scenario of IGW as brought out by 
the temperature dependence of z/. 

Let £g = Pg^o denote the dimensionless energy per 
nbNN contact at the growth tepmperature Pq^. Then, 
an A^-stcp IGW configuration, C, having a total of Nc{C) 
such contacts will have an energy, Eg{C) = £gNc{C). As 
illustrated in Fig. 1(b) and 1(c), configurations with the 
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same energy are generated with different probabilities. 
We may rewrite the growth probabihty, Pigw{^,C), as 
follows. 

N 

PiGw{N;C) = [|pj(r,-;ro,ri,...,r,_i) (2) 

where PsAw{N) = z~^{z — l)~'^'^~'^^ is the probability of 
generating an A''-step SAW configuration and £{C) is the 
energy per contact to be assigned to the configuration 
if it were to be considered as a member of a canonical 
ensemble. 

1 ^ 

^^^^ = ~ lfe(r,;ro,ri, ...,r,-_i)] (4) 

It is now clear that different configurations with the same 
number contacts could be assigned different values of 
£{C) because their growth probabilities are different. In 
other words, for a given value of the growth parameter, 
£g, the mapping of IGW to IS AW gives rise to a distri- 
bution of the dimensionless energy per contact, £. 

Assuming that eo is a constant, a distribution in £ 
corresponds to a distribution in /3. This implies that the 
IGW configurations grown at a given temperature 
can be considered as ISAW configurations, but sampled 
at temperatures drawn from a distribution in /3. We have 
discussed this recently for IGW on a honeycomb lattice 
[8]. We have shown that a sharply peaked distribution in 
/3 can be associated with any given /3g > (the broadest 
distribution, numerically obtained for Pg = oo, peaks at 
f3 ~ 1.21 with a FWHM ~ 0.03). In the athermal limit 
{(3g = 0), the IGW corresponds to ISAW at a unique 
temperature given by /3 = log 2, a result obtained first by 
Poole et al [9]. Since the distribution in /3 is sharp, the 
peak value may be taken to provide a well defined canon- 
ical or 'bath' temperature at which most of the IGW 
configurations can be considered as ISAW configurations. 
The ones that correspond to different temperatures will 
have to be equilibrated at the peak temperature. 

Alternatively, if IGW were to be considered as an 
ISAW, then it should represent an equilibrium config- 
uration at a uniquely defined 'bath' temperature. We fix 
the bath temperature, j3, by assuming that the peak po- 
sition of the distribution in £ can be identified with /3eo- 
There is no a priori reason to assiimc that the average en- 
ergy per contact for the equilibrium configuration should 
be the same as eo, a parameter introduced for sampling 
the locally available sites during its growth. Hence, the 
distribution in £ can be taken to be proportional to a 
distribution in e, peaking at eo- 
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FIG. 4. Specific heat as a function of bath temperature, 
T = f3~^ . The sharp peak at T ~ 1 corresponds to Pa ~ 4.5, 
and hence to the ^-collapse transition. The continuous hnc is 
a guide to the eye. Inset: Inverse of bath temperature, /3, as 
a function of the inverse of growth temperature, Pa- 

We have obtained the bath temperature, P{N), and 
the width, <j(N), of the distribution in e as a function of 
N for a given Pg , basically from the first and second mo- 
ments of the distribution in £. Then, we have estimated 
their asymptotic values by fitting them to a simple form, 
y(N) = y + {A/ N^) where y(= P or cr), A and B are 
adjustable parameters. We have presented the estimated 
P values as a function of Pg in the inset of Fig. 4. We 
find that the full range of Pg G [0, oo] is mapped into a 
narrow range of bath temperatures, P € 0.42, ^ 1.12] 
(e [log 2,^ 1.2], on honeycomb lattice [8]). It may be 
noted that the 6'-point, Pg ^ 4.5, corresponds to /3 ^ 1. 

From the asymptotic variances, (t^(/3), wc have ob- 
tained the specific heat per contact, c(/3) = P'^a'^iP), and 
presented them in Fig. 4 as a function of the bath tem- 
perature P^^. The sharp peak seen at about /? ~ 1 corre- 
sponds to the collapse transition at the 6'-point. This, in 
fact, validates the view that a definite bath temperature 
can be associated with the IGW. 

But, there is no known transition that can be associ- 
ated with the excess specific heat seen as a broad hump 
above the ^^-peak, because this region is in the SAW phase 
as far as the universal exponents are concerned (Fig. 3). It 
is therefore of interest to understand what is responsible 
for this excess specific heat. Recently, hyperquenched 
glasses have been shown [10] to exhibit excess specific 
heat (Fig. 4 of Ref.[10]), strikingly similar to what we 
have observed for the IGW (Fig. 4) above the 0-point. 
The dimensionless energy per contact, £{C), defined in 
Eqn.4, is indeed an average of such values that can be 
evaluated during the growth process. This implies that 
a distribution of £ can be associated with every configu- 
ration generated. Moreover, the IGW configurations are 
clearly much more compact (see Fig.l of ref.[6]) than the 
typical SAWs belonging to the same universality class. 
It is therefore reasonable to consider them as "frozen" 
globules. 
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FIG. 5. A schematic illustration of how the growth of an 
IGW can be viewed as a hierarchical process. The configura- 
tions are coded as strings of Os, Is, 2s and 3s, enclosed within 
square brackets, where the labels 0, 1,2 and 3 correspond to 
steps in the +x, — y, —x and +y directions respectively. The 
various paths in the hierarchy are taken with different proba- 
bilities (see text). The tree is constructed in such a way that 
the final configurations are numbered in increasing order from 
left to right. Shown just below the tree is the energy land- 
scape for all the 5-step walks whose first step is along the +x 
direction. Of course, the probability of realising a point on 
the landscape depends on the growth temperature, /3q^. The 
global minimum energy {— —2) configurations are indicated 
by their respective codes. And below this is a schematic pic- 
ture of the energy landscape for asymptotically long walks. In 
the case of IGW, the number of available (or realisable) final 
configurations decreases as the walk proceeds to grow. This is 
illustrated by shaded regions becoming progressively darker. 
Exactly which point on the landscape is finally reached is de- 
cided by the value of /3q^. In the case of ISAW, however, 
all the configurations having their energies within an interval 
(schematically indicated by the shaded region) determined by 
f3~^ will be sampled. 

It may be noted that the correspondence between /?g 
and /? whose existence is dictated by Eqn.(4) forms the 
basis of this study. And, the fact that the full range 
of /3g S [0,oo] maps into a finite range of canonical 
/? G ^ [0.42,1.12] has subtle physical implications. For 
example, as depicted in Fig. 5, the growth of an IGW 
can also be considered as a hierarchical approach to- 
wards realizing a particular configuration. Every step 
taken reduces the number of available configurations, or 
equivalently, restricts the accessible region of the energy 
landscape in a progressive manner. This implies that ir- 
reversible growth is equivalent to breaking the ergodicity 
of the system. The probability of taking a certain path 
in the hierarchy depends on the timing parameter, (3g. 
On the other hand, in the canonical ensemble picture, we 
sample all the configurations whose energies lie within an 
interval defined by the bath temperature, /3~^ (schemat- 



ically illustrated in Fig. 5). In particular, we expect to 
sample only those configurations with global minimum 
energy when = 0. In contrast, with = 0, the 
IGW algorithm will generate a few zero energy (ather- 
mal) configurations as well, besides those with global 
minimum energy; hence, the corresponding /3~^ will be 
greater than zero. And, larger the value of the coordina- 
tion number, z, of the lattice, smaller will be the num- 
ber of such athermal configurations and hence larger will 
be the value of l3~'^ to which it corresponds. Similarly, 
the distribution of NN contacts for the IGW configura- 
tions generated at /3g = deviates from that obtainable 
for SAW, and hence the corresponding /3 will be an z- 
dependent nonzero value. 

In summary, we have shown that the IGW configura- 
tions can be considered as members of a canonical en- 
semble {i.e., as ISAW configurations) if the energy per 
contact can be considered as a random variable. In 
general, a meaningful statistical mechanical description 
of an irreversible growth process involves an element of 
self-generated disorder. The signature of this is seen as 
a broad hump in the specific heat above the 6- point. 
That these configurations are generated in an hierarchi- 
cal manner, as implied by the specific growth rule, pro- 
vides additional support to the conjecture that they may 
be taken to represent hyperquenched polymer configu- 
rations. Conformational dynamics of IGW could throw 
further light on this conjecture. In fact, the IGW seems 
to illustrate the generic possibility of a growth process 
giving rise to hyperquenched states of a system, if it is 
faster than the configurational relaxation. 
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